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ABSTRACT 

We present the first 7-millimetre observations of two high-redshift, Lya-bright radio 
galaxies (MRC 2104-242 and MRC 0943-242) performed with the 2 x 2 GHz instanta- 
neous bandwidth of the Compact Array Broadband Backend (CABB) at the Australia 
Telescope Compact Array (ATCA). The aim was to search for CO (1-0) emission in 
these systems and test the millimetre capabilities of CABB for performing spectral 
line observations at high redshifts. We show that the stable band and enhanced veloc- 
ity coverage of CABB, combined with hybrid array configurations, provide the ATCA 
with excellent 7-mm capabilities that allow reliable searches for the ground transition 
of CO at high redshifts. In this paper we explicitly discuss the calibration procedures 
used to reach our results. We set a firm upper limit to the mass of molecular gas in 
MRC 2104-242 [z = 2.5) of Mh2 < 2 x 10^° (ax/0.8) Mg. For MRC 0943-242 (z = 2.9) 
we derive an upper limit of Mh2 < 6 x (q;x/0.8)Mq. We also find a tentative 
3(j CO detection in the outer part of the giant Lya halo that surrounds MRC 0943- 
242. The 30-33 GHz radio continuum of MRC 2104-242 and MRC 0943-242 is rehably 
detected. Both radio sources show a spectral index of a « —1.5 between 1.4 and 30 
GHz, with no evidence for spectral curvature within this range of frequencies. 

Key words: galaxies: high-redshift - galaxies: active - galaxies: ISM - galaxies: 
individual: MRC 2104-242 - galaxies: individual: MRC 0943-242 - techniques: inter- 
ferometric 



1 INTRODUCTION 

Cold gas is a primary component in galaxy formation pro- 
cesses such as star formation and disk growth. However, de- 
spite detailed studies of cold gas in the nearby Universe, it 
is still difficult to trace similar qu antities of cold gas be yond 
our Galactic backyard. Recently, iTacconi et al.l (|2010l l and 
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iDaddi et all l|2010l ) observed that star-forming galaxies at 
high redshifts are likely to contain a much larger fraction of 
their total mass in the form of molecular gas compared with 
nearby massive spiral galaxies. Recent simulations support 
this idea that the molecular gas co ntent of galaxies increases 
when g oing to higher redshifts ijObreschkow fc Rawlingd 
l2009al lbl: lObreschkow et al.1 l2009al lbir These results demon- 
strate that extensive studies of cold molecular gas in the 
early Universe are becoming feasible with existing radio tele- 
scopes. 
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Powerful radio galaxies enable comprehensive studies 
of the cold ISM throughout the Universe. Their strong 
radio sources provide a background continuum against 
which we can search fo r foreground neut r al and molecular 
gas i n absorption (e.g. Uson et al. 199ll : Vermeulen et al.l 



I2OO3I : iKanekar et al.1 EoOTI : I Carilli et al while their 

host galaxies are generally in a very specific stage of 
galaxy evolution. Detailed studies at low and interme- 
diate redshifts reveal that powerful radio galaxies are 
frequently asso c iated with gas-rich galax y mergers (e.g. 
iHeckman erahl Il986l : [b aum et ahl Il992|'). often contain 
youn g stellar populati o ns ([Tadhunter et al]|2005l : iHolt et al] 



I2OO7I : iLabiano et al] l2008h. an d ma ny display strong 



jet-ISM interactions l|Tadhunteii 11 991: Vi llar-Martm et al 
19991: [Clark et al.lll998l:lEmonts et al.„2005nMorganti et al 



2005al lbl: iHoh et al.l l2008h . At high redshifts {z > 2), lu 



minous radio galaxies (Lsoomhz > 10^^ W Hz ^) are 
among the most massive galaxies in the early Universe (see 
iMilev fc De Breucklliooi . for a review). They are typically 
surrounded by proto-clusters, whic h are thought to be the 
ancestors of rich local clusters (e.g. iPentericci et al.ll2000al : 
IVenemans et alj|2007l ). The high-z radio galaxies and sur- 
rounding proto-cl uster gas and galaxies often interact with 
one another (e.g. iNesvadba et al.l l2008l : llvison et al] I2OO8I ') 
and are therefore laboratories for studying the formation 
and evolution of galaxies and clusters as well as investigat- 
ing the relationship between early star formation and AGN 
activity. 



Since iBrown fc Vanden Bmitl |l993) first observed CO 
gas (the strongest tracer for molecular hydrogen) at a red- 
shift beyond z = 2, intensive searches for CO in high- 
z radio ga laxies d uring the early 1990s were unsuccessful 
vans et al. 1996. : Ivan Oiik et al. 1997). Since then, stud- 
ies of individual radio galeixies at 2; ~ 2 — 5 with synthesis 
radio telescopes have found CO emission (tracing molecu- 
lar gas masses of a few xlO^° - lO^Mp)) in a numb er of 



these systems (e.g. Scov illc ct al. 
2OOOI: iDe Breuck et al.l [2003a . b, 



Nesvadba et 



„ ^ . 

I 2005h : lOmon? (|2007l ): iMilev fc De Breuckl ||2008D for 



200S 



1997; 
2005: 



see also 



Papadopoulos ct al] 
Klamcr ct al. 20051: 



Solomon fc Vanden Bout! 



views). In some cases CO is observed to be resolved on scales 
of several tens of kpc (e.g. IPapadopoulos et al.ll200ol l. This 
indicates that large amounts of cold molecular gas may be 
relatively common in high-z radio galaxies. However, the 
major observational limitations for starting comprehensive 
studies of CO in high-z radio galaxies have been the very 
limited velocity coverage of existing mm-spectrometers (of- 
ten not much wider than the velocity range of the CO gas 
and/or the accuracy of the redshift) plus the fact that most 
observatories can only target the higher order rotational 
transitions of ^^CO. 



Although the higher order CO lines are likely to have 
a higher flux density than the lower ones in the nuclear 
starburs t/AGN regions, whe r e gas is den se and thermally 
excited, IPapadopoulos et al. I (|2000l . l200lf ) suggest that the 
opposite may be true for large reservoirs of less dense and 
sub-thermally excited gas that is more widely distributed. In 
fact, various studies of the low-order CO transitions in differ- 
ent types of high-jz galaxies reveal molecular gas that is sub- 



thermally excite tfl (Grov e et al]|2003l: iHainline et al]|2006l : 
iDannerbauer et al.„200a:iRiechers et al.ll2010f) or distributed 
in extended reservoirs jPaddi et al.ll2010l : ICarilh et al.ll2010l : 
llvison et al] |2010| . |2011^ . Cold CO gas distributed across 
the host galaxy may thus be much easier to detect in the 
lower CO transitions than generally assumed from studies 
of the higher transitions. Moreover, with uncertainties in ex- 
citation properties of the gas, observations of the rotational 
ground-transition of the CO molecule [^^CO(l-O) - referred 
to as CO in the remainder of this paper] provide the most 
accurate mass estimate of the overall molecular gas content 
in these systems. 

Since April 2009, the Australia Telescope Compact Ar- 
ray (ATCA) has a new broad-band backend system (the 
Compact Array Broadband Backend or CABB). CABB of- 
fers an instantaneous bandwidth of 4 GHz, split over 2x2 
GHz observing bands, both with all Stokes polarisation 
parameters a nd 2048 channe l s (i.e. spectral resolution of 
1 MHz); see (Ferris fc WilsonI l|2002l ): Iwilson'et al.1 (|201ll) . 
ATCA/CABB has millimetre observing capabilities at 3mm 
(83.9 - 104.8 GHz), 7mm (30.0 - 50.0 GHz) and 15mm 
(16.0 — 25.0 GHz). This, in combination with hybrid array 
configurations with baselines as short as 31m, makes the 
upgraded ATCA an excellent facility to detect and spatially 
resolve molecular gas in high-z radio galaxies by targeting 
the lower rotational CO transitions (see Sect.[2]for more de- 
tails) . A remarkable example of this is the recent detection 
of C0(2-l) in the distant (z = 4.8) sub - millim etre galaxy 
LESS J033229.4-275619 bv ICoppin erahl l|2010t ). 

To test the spectral-line performance of CABB over the 
2x2 GHz bandwidth, we used the 7mm band to search 
for CO(l-O) in two high-z radio galaxies from th e Molon- 
glo Reference Catalogue (|McCarthv et al] Il990l ). namely 
MRC 2104-242 (2 = 2.5) and MRC 0943-242 {z = 2.9). 
These two sources are part of a larger sample of high-z ra- 
dio galaxies that we aim to target with CABB in order to 
perform a systematic search for CO(l-O) in these systems. 

MRC 2104-242 and MRC 0943-242 both have a redshift 
that corresponds to a critical epoch in galaxy formation {z ~ 
2.5 — 3), at w hich there is a dramatic increase in sub-m m flux 
ll Archibald et al. 2001; Small ct al. 2002; Chapma n et al.1 
I2OO5I ) and the s pace-dens i ty of (radio-loud) quasars reaches a 



s pa ,_-^-r 

maximum (e.g. |Pei|[l995l : [ Shaver et al.ll 19961: iRichards et al.l 
I2OO6I ). HST observations bv IPentericci et all (120011 ) show 
that MRC 2104-242 and MRC 0943-242 both have an optical 
continuum that i s clumpy and elongated in the direction o f 
the radio source l|Pentericci et al.ll2000bl : ICariUi et al]|l997l ) . 
IViUar -Martin et al.l (|2003l ) show that they both contain a 
giant Lya-halo (^ 100 kpc in diameter). For MRC 2104- 
242 the Lya gas is distributed roughly along the radio axis 
in what appear s to be a rotating structu re with a diame- 
ter > 120 kpc (|Villar-Martm et aLllioO^ . MRC 0943-242 
shows a quiesce nt Lya-halo that extends well beyond the 
radio structure l|ViUar-Martm et al.ll2003l). MRC 0943-242 
also shows a deep Lya absorption l|Rottgering et al.|[T995l : 
IJarvis et al.l |2003| ) , indicating that large amounts of neu- 



^ [Harris et al. I 1 I2OIOI) and iDanielson et aL I | j2010l) point out that 
a multi-component inter-stellar medium, rather than sub-thermal 
excitation, may better reflect the physical properties of the molec- 
ular gas in high-2 systems. 
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Table 1. Observations 



-CO(I-O) 
-C0(2-1) 
CO(3-2) 
-C0(4-3) 
-CO(5-4) 
-CO(6-5) 



Redshift 

Figure 1. Observing frequency of the first six rotational tran- 
sitions of CO that can be targeted with the CABB millimetre 
system (3, 7 and 15mm) plotted against the redshift of the CO- 
emitters. 



Source 


Array 


Obs. date 


tint (h) t'cciitial (GH 


MRC 2104-242 


H75 


09JUL21 


6.99 


33.000 






09JUL23 


4.15 


33.020 




H168 


09SEP18 


3.00 


33.000 






09SEP19 


3.14 


33.000 






09SEP20 


2.19 


33.000 


MRC 0943-242 


H75 


09JUL11 


3.10 


30.001 






09JUL12 


3.20 


30.010 






09JUL14 


2.88 


30.001 




H168 


09MAY10 


3.29 


30.001 






09MAY11 


2.68 


30.015 






09MAY12 


3.18 


30.001 


Notes - tint is the effective 


on-source integration time ( 


.0. not 



including overheads). 



tral gas are present in this system. From fitting tiie spec- 
tral energy dis t ributi ons of the host galaxies with Spitzer, 
ISevmour et"ai] (|2007l ') derive a total stellar mass of a few 
xlO^^M0 for both systems. 

In Sect. [2] we present our CO observations and explain 
in more detail the enhanced capabilities of the ATCA for 
studying molecular gas at high redshifts. Section |3] shows the 
result regarding both the performance of CABB for doing 
these high-z CO studies as well as the scientific outcome of 
our observations of MRC 2104-242 and MRC 0943-242. In 
Sect. |4] and [5] we discuss the scientific results and conclude 
that the upgraded ATCA is a world-class facility for spectral 
line observations of the cold molecular gas at high-2. 



2 OBSERVATIONS 

During the period May - September 2009, MRC 2104-242 
and MRC 0943-242 were observed with ATCA/CABB. De- 
tails of the observations are given in Table [1] 

Figure [T] shows the observing windows for the vari- 
ous transitions of extra-galactic CO currently available with 
CABB in the 3, 7 and 15 milhmetre bands. For MRC 2104- 
242 and MRC 0943-242 we targeted the ground-transition 
CO(l-O) with the CABB 7mm system. The redshift of 
MRC 2104-242 {z = 2.491) corresponds to an observing 
frequency of 33.0 GHz, which is also one of the optimum 
CABB frequencies for centringthe 7mm bandQ The redshift 
of MRC 0943-242 (z = 2. 9185l3 corresponds to an observing 
frequency of 29.4 GHz, which is outside the nominal 7mm 
CABB band (30.0-50.0 GHz). Nevertheless, when centring 
the band at 30.001 GHz, data is obtained down to 29 GHz. 
Observations of MRC 0943-242 therefore served as a good 
test of how well the CABB system behaves at the very edge 
of the 7mm band. 

The total observing time for each source - including 
overhead and calibration - was roughly 40 hours (see Table 
[T]and Sect. 12. 1| ). The observations were spread over the two 

^ See the online ATCA Users Guide for details: |http: //www[ 
narrabri.atnf.csiro.au/observing/users_guide/html/atug.html 
^ This redshift corresponds to the velocity of the mo st prominent 
H I a bsorption in the Lya profile of MRC 0943-242 l lJarvis et al.l 
I2OO3I) : see Sect. [2X3] for more details. 



most compact hybrid array configurations (H75 and H168) 
in order to minimise the effect of atmospheric phase fluctua- 
tions (which worsen with increase in baseline length: iKlameJ 
l2006t ). Both array configurations include five antennas that 
are spread across both an east-west as well as a north-south 
spur. This ensures a reasonably good w^^coverage (see Fig[2} 
even from observing runs as short as 6-8h, during which a 
source is targeted only above an elevation of 30° (in order 
to avoid high airmasses, which d ramatically increase the 
system temperature; see iKlameij I2OO6 ) . The system tem- 
peratures ranged from 70 — 132K for the observations of 
MRC 2104-242 and from 130 - 180K for MRC 0943-2420 

At the high central frequencies of our observations 
(30-33 GHz), the coarsest spectral-line mode of CABB 
(2x2 GHz bands with 1 MHz spectral resolution) provides 
a velocity resolution of 9.5-11 km s~^ across an effective ve- 
locity coverage of at least 17,000 km s^^ per 2 GHz band. 
Both 2 GHz bands were centred around the same observing 
frequency, but because they are not mutually independent, 
only one 2 GHz band was used in the final data analysis. 

Observations were done as much as possible during the 
night and in good weather conditions (to avoid decorrelation 
due to atmospheric phase instabilities) and only above 30° 
elevation (to avoid high system temperatures due to large 
airmasses). Because our observations also served as a test 
for the performance of CABB for spectral-line observations, 
in the following Section we will explain the details of several 
crucial calibration steps. For the data reduction and analysis 
we used MIRIAD and KARMA. 



2.1 Calibration, overheads and data reduction 

Our general observing strategy was as follows: a strong cal- 
ibrator was observed at least three times during each run 
in order to check the reliability of the bandpass calibration. 
A secondary (phase/gain) calibrator was observed roughly 
every 10 minutes. Flux calibration was done at least once 
during each run. Pointing solutions of the antennas were 
checked and updated every hour, or every time the telescope 



* More details on theoretical estimates of Tgys values at 7mm 
can be found in the online ATCA Users Guide. 
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MRC 2104-242 (33 GHz) MRC 0943-242 (30 GHz) 




u Cki) u (lit) 



Figure 2. The lii^coverage of the combined observing runs of both MRC 2104-242 and MRC 0943-242 (Table[TJ, using the hybrid H75 
and H168 arrays. 



slewed more than ~20° on the sky. Taking into considera- 
tion the conservative nature of this cahbration strategy, the 
overheads due to cahbration and slewing were about 50% 

2.1.1 Phase/gain calibration 

For phase calibration we performed a 2 minute scan on a cal- 
ibrator close to our target source roughly every 10 minutes, 
although target scans were decreased to 5 minutes in poor 
weather conditions and increased to 15 minutes when atmo- 
spheric phase stability was excellent. For MRC 2104-242 we 
used PKS B2008-159, PKS B2128-123 or PKS 2149-306 as 
phase calibrator. For MRC 0943-242 we used PKS 0919-260. 
Phase calibration was done in a standard way. 

2.1.2 Bandpass calibration 

In order to test the quality of the bandpass calibration at 
7mm across the full 2 GHz band, we observed a strong cal- 
ibrator (PKS B0537-441, PKS B1253-055, PKS B1334-127, 
PKS B1921-293 or PKS B2223-052) at least three times dur- 
ing each run (unless the run was cut short due to weather). 
We noticed that weather and atmospheric conditions at the 
ATCA site can introduce frequency dependent temporal 
gain fluctuations across the wide CABB band, which can 
have a significant effect on the quality of the bandpass cal- 
ibration at 7mm. It is therefore essential to obtain at least 
one good scan on the bandpass calibrator during good at- 
mospheric conditions. For MRC 0943-242 we chose the best 
quality bandpass calibrator scan for calibrating our data. In 
case more than one bandpass calibrator scan was deemed 
suitable, we applied the bandpass solutions to that part of 
the data observed closest in time to the respective calibrator. 

For MRC 2104-242 the strong phase calibrators 
PKS B2008-159 and PKS B2128-123 (with observed fluxes 

^ We estimate that in order to reach the potential maximum 
efficiency with less conservative cahbration, overheads should be 
considered to be at least 30%. 



of F33 GHz ~ 1.9 and 1.8 Jy respectively) were suitable for 
bandpass calibration. This allowed us to obtain a bandpass 
solution roughly every 10-15 minutes. We used a new feature 
in the MIRIAD task mfcal to interpolate between consec- 
utive bandpass solutions in order to compensate for possible 
frequency dependent gain variations that slowly fluctuate in 
time. 



2.1.3 Flux calibration 

For MRC 2104-242, flux calibration was done by observing 
Uranus at the time that it was at roughly the same eleva- 
tion as the phase calibrator and target source during each 
run. The presence of a weak radio continuum from the lobe- 
dominated high-z radio galaxies in our 7mm data (which 
are not expected to significantly change their fiux densities 
over time-scales of a few months) allowed us to compare the 
relative flux calibration between the various runs, which re- 
mained constant within 15%. Our absolute flux calibration 
used the available MIRIAD-model for Uranus. This model 
did not take into account changes in the planet's orientation, 
which introduce time-v ariations of up to 10% in its bright- 
ness temperature (see iKramer et al ] l2008l : IWeiland et al.l 
l2010t ). potentially leading to a significant error in absolute 
flux calibration. During one of the runs we also observed 
PKS B1934-638, which confirmed our Uranus-based abso- 
lute flux calibration to an accuracy of ~18%. We therefore 
estimate the overall (relative -I- absolute) uncertainty in the 
flux calibration of MRC 2104-242 to be within 30%. 

For MRC 0943-242, Uranus was not visible during our 
observing runs. For flux calibration we therefore observed 
the uhra-comp ac t HII region G309 [G309. 9206-1-00.4790; 
lUrauhart et~al] |20o3), with our pointing centred at 
RA(J2000) = 13:50:42.35, dec(J2000)=-61:35:09.78] when it 
was at roughly the same elevation as the phase calibrator. 
We calibrated the flux of G309 against Uranus, which we ob- 
served roughly half a day later for each run. The flux of G309 
was stable over our six observing epochs and the relative flux 
calibration between the six different runs was within 13%. 
From our data we derive a value of Sao ghz ~ 1.31 ± 0.07 Jy 
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for the sh ortest baselines at wh ich the source is unresolved. 
Recently, iMurphv et all (|2010l ') derived a flux density of 
S32 GHz = l.liO.llJy for G309, also using Uranus as 
flux calibrator. In order to verify the accuracy of our ab- 
solute flux calibration, we observed PKS BI934-638 during 
three of our observing epochs. When using PKS B 1934-638 
as flux calibrator instead of Uranus, the absolute fluxes de- 
rived from our data are on average ~ 15 % lower. This un- 
certainty in absolute flux calibration is consistent with the 
difference between our flux estimate for G 309 (which we 
used t o calibrate our data) and that made bv lMurphv et all 
This may again reflects variations in the brightness 
of Uranus that were not accounted for by the existing mod- 
els (see previous paragraph). The spectral index of G309 
changes at most a few per cent across th e 2 GH z band at 
30 GIfz, in agreement with iMurphv et al] (|201Cf ). In all, we 
therefore estimate that for MRC 0943-242 the overall (rela- 
tive -I- absolute) uncertainty in our flux calibration is within 
30%. 

After flagging and bandpass, gain and flux calibration, we 
subtracted the continuum from the line data in the uv- 
domain by applying a linear fit to the channels across the 
full 2 GHz band (for MRC 0943-242 we excluded from this 
fit the channels in which we found a tentative CO signal, see 
Sect. 13.31 although this has no significant effect when fitting 
the full 2 GH z band). Subsequently, a robust -|-1 weighted 
jBriggd Il995l ) continuum map and line data set were cre- 
ated by Fourier transforming the Mi>-data and, in case of the 
continuum map, cleaning the signal. We then translated the 
velocity axis to match the optical, barycentric rest-frame 
velocity at the redshift of MRC 2104-242 and MRC 0943- 
242. The red shift of MRC 2104-2 42 (z = 2.491) has been 
confir med bv lOverzier et ahl |20o3) and I Villar-Martfn et al.l 
l|2003h through observations of the Lya and various metal 
emission-lines. For MRC 0943-242, we chose to centre our 
observations at the redshift of th e prominent H I ab sorption 
in the Lya profile {z = 2.9185; Ijarvis et al]l2003l '). which 
likely represe nts the bulk of the cold neutral gas in this sys- 
tem (see also lDe Breuck et al.ll2003lj ). Table [2] shows details 
of the final data products that we obtained from our obser- 
vations (some of these are described further in Sect. [3}. 



3 RESULTS 

3.1 CABB performance 

Figure O shows the 33 GHz radio continuum map of 
MRC 2104-242 and a spectrum at the location of the centre 
of the host galaxy. The continuum image has an rms noise 
level of 29 /xjy beam~^ (after tint ~ 19. 5h; see Table (2]), 
demonstrating the effectiveness of ATCA/CABB for deep 
millimetre continuum studies. The 2 GHz spectrum has a 
large velocity coverage of ~ 17,000 km s~^ with an rms 
noise m each 1 MHz channel (Av = 9.6 km s"^) of cr = 0.45 
mjy bm~^, with no significant systematic bandpass effects. 

Figure [J] shows the 30 GHz radio continuum map of 
MRC 0943-242 (with an rms noise of 33 ^iJy beam"^) and 
an off-nuclear spectral line profile. In this case, at the edge 
of the 7mm band, half the observing band lies outside the 
nominal CABB range (Sect.[2|, where there are instrumen- 
tal low-level structures in the noise or in the bandpass at 



Table 2. Data 









Effective int. time (h) 


19.5 


18.3 


Target frequency (GHz) 


33.02 


29.417 


Redshift 


2.491 


2.9185 


Beam size (arcsecxarcsec) 


11.7 X 7.6 


11.5 X 9.0 


Beam PA [PA (°)] 


97.2 


87.5 


Av (km s"-"-) 


9.6 


11.0 


CTcont (AiJy bm"^) 


29 


33 


5cont (mJy) 


4.0 


3.3 


'''line (mJy bm~^ch~^) 


0.45 


0.90 


L^Q (K km s-lpc2) 


< 2.6 X IQi'^ 


< 7.3 X 



Notes - Effective int. time is the total effective on-source 
integration time of all runs in both configurations combined. 
Target frequency (GHz) is the observing frequency of the 
expected CO(l-O) line at the redshift of our sources (see text for 
details). Av (km s~^) is the velocity resolution per 1 MHz 
channel. iTcont is the rms noise level of the continuum image 
after tint ■ Scont is the integrated continuum flux of the radio 
source at the target frequency, crnne is the rms noise level of the 
full-resolution line data per 1 MHz channel after tint- L'qq gives 
the upper limit on the CO luminosity (see text for details). 

about the Icr level of the full (IMHz) resolution data (Fig. 
U)). In addition, the noise starts to vary beyond the nomi- 
nal CABB range (Fig. 14]:). After an effective on-source in- 
tegration time of 18. 3h, we derive a noise level at 29.4 GHz 
of cr = 0.9 mJy bm~^ per 1 MHz channel (Av = 11 km 
s~^), i.e. twice the noise level at the optimum observing fre- 
quency of 33 GHz (see above). However, as can be seen in 
the Hanning-smoothed data of Fig. 2]:, the noise level peaks 
at our target frequency of 29.4 GHz and is significantly lower 
throughout most part of the band, even below the nominal 
edge of 30 GHz. We therefore conclude that up to ~ 0.8 
GHz below the nominal 7mm band, CABB is still suitable 
for spectral-line work. 

ICoppin et all (|2010D detected C0(2-l) in a 2: = 4.8 sub- 
millimetre galaxy, which was observed with CABB at 40.0 
GHz (i.e. towards the other end of the 7mm band compared 
to our 33/30 GHz observations). They find noise levels of 
a « 0.44 mJy bm~^ per 1 MHz channel and a bandpass 
stable enough to detect their CO signal at about the 5(t 
level when binning across >10 channels. The data quality 
at 40 GHz thus appears comparable to that at 33 GHz as 
presented in this paper, giving a good indication for the ex- 
cellent performance of CABB across the entire ATCA 7mm 
band. 

3.2 MRC 2104-242 

MRC 2104-242 (Fig. Eli) is resolved at 33 GHz with a to- 
tal flux of 3.5 mJy. The continuum structure consists of 
two components on either side of the optical host galaxy, 
in agreement with 4.7 and 8.2 GHz VLA continuum ob- 
servations that identifie d it as a double lobed radio source 
ijPentericci et all l2000bl ). The bright northern lobe has a 
peak flux density of iSaaGHz = 2.7 mJy bm~^, while the 
fainter southern lobe has Sssghz ~ 0.39 mJy bm^^. Even at 
33 GHz the radio continuum structure is dominated by the 
radio lobes and no core component (at the location of the 
optical nucleus) is seen in our data. We set a conservative 
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upper limit to the 33 GHz core flux density of S'core-ssGHz < 
0.4 mjy bm~^. Figure [5] shows that the integrated flux of 
MRC 2104-242 has a steep spectral index between 1.4 GHz 
and 33.0 GHz, with a = —1.56 (where Fv oc There is no 
evidence for spectral curvature within this range of frequen- 
cies. This is in agreement with spectral index o bservations of 
high-g ultra-steep spectrum radio sources bv iKlamer et al.l 
(|2006l l. 

No CO is detected in MRC 2104-242, either at the loca- 
tion of the host galaxy or at the position of the radio source. 
We derive a firm upper limit on the CO emission-line lumi- 
nosity in MRC 2104-242 by assuming a potential 3(t signal 

S using 



ScoAV = 3(7 Av 



500 km s- 



Jy ■ km s 



(1) 



with a the noise level per 1 MHz channel in one beam (in Jy) 
and Av the width of one 1 MHz channel (in km s~^). The CO 
luminosity (upper limit) can t hen be calculated following 
ISolomon fc Vanden BoutI (|2005l . and references therein): 



2 / i^rcst . 



GHz' 



^(1 + ^) 



(2) 



with L'qq expressed in K km/s pc^ and with _Dl ~ 20018 
Mpc th e luminosity distance of MRC 2104-242 (following 
IWrightlliooa fl. For MRC 2104-242, ScoAV < 0.094 Jy km 
s"\ hence L'co < 2.6 x lO"' K km s"^ pc^. 



10= 



10' 
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a = -1.56 



10° V- 

10° 
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10' 

Observed Fiequency (GHz) 



10= 



10= 



E 



10' 



MRC 0943-242 

a = -1.44 



3.3 MRC 0943-242 

The radio source MRC 0943-242 (Fig. [Dj) has a flux den- 
sity of 3.3 mJy bm~^ and is unresolved in our data. Higher 
resolution continuu m observations at 4.7 and 8.2 GHz by 
ICariUi et al.l (| 19971 ) show that the radio source consists of 
two lobes that are separated by 4 arcsec. When compar- 
ing the flux of our 30 GHz data with the integrated flux at 
1.5, 4.7 and 8.2 GHz (jCarilh et al.lll997l ). Figure [5] shows 
that MRC 0943-242 has a steep spectral index between 1.5 
GHz and 30 GHz with a = —1.44. Similar to the case of 
MRC 2104-242, there is no evidence for spectral curvature 
within this range of frequencies. 

No CO is detected at the central (nuclear) location of 
MRC 0943-242. When estimating an upper limit on L'qq 
in MRC 0943-242 (potential 3ct detection smoothed across 
500 km s"^), we derive L'co < 7.3 x 10^" K km s"^ pc^ (for 
Dl = 24242 Mpc, which corresponds to a angular- size scale 
of 7.65 kpc/arsec for MRC 0943-242: IWrightlf2006l ). 

3.3.1 Tentative off-nuclear CO detection 

As can be seen in Fig.[3]b,d,e, we find a tentative, off-nuclear 
3(T CO (1-0) detection in the Hanning smoothed data of 
MRC 0943-242 (with a the noise level at the frequency that 
corresponds to the tentative detection, see the arrow in Fig. 
3]:). The tentative CO signal spreads over an area about the 



^ See I http://www.astro.ucla.edu/~wright/CosmoCalc.htmll for 
Ned Wright's online cosmology calculator that we used to deriv- 
ing luminosity and angular-size distances. Throughout this paper 
we use Ho = 71 kms~^ Mpc^^, f2M = 0.3 and Ha = 0.7. 



10" 10' 10= 

Observed Frequency (GHz) 

Figure 5. Spectral index of MRC 2104-242 (top) and MRC 0943- 
242 (bottom). Shown is the total i ntegrated flu x . Da ta at 
1.4, 4.7 and 8.2 GHz are taken from ICarilli et al.l lll997l ') and 
iPentericci et al. ] I^OOOb), with errors corres ponding to 2% of the 
flux at these wavelengths (as estimated bv lCarilli et al.lll997t) . 

size of one synthesised beam roughly 60 kpc NE of the cen- 
tre of the host galaxy. It peaks at v ^ —100 km s~^ with 
a flux density of ~ 1.8 mJy bm~^(with a tentative second 
peak present at the 1 mJy bm~^ level around v ~ —500 
km s~^). The estimated luminosity of the tentative double- 
peaked CO signal is L'qq ~ 8 x 10^° K km s~^ pc^ (Equa- 
tion[2)|. Both the H75 and H168 array data show indications 
for this tentative CO signal (Fig. |4^). However, because of 
the low-level (la) structure in the noise/bandpass beyond 
the nominal 7mm observing band (see Sect. 13. l| . our results 
did not improve by further smoothing/binning the data in 
velocity. Our tentative 3cr detection thus needs to be veri- 
fied with additional observations before conclusions can be 
drawn. 

3.4 The environments of high-z radio galaxies 

The large instantaneous velocity coverage of CABB (see 
Sect. [21 also makes it possible to search for CO emitters 
in the field of our high-z radio galaxies. The full width half 
maximum (FWHM) of the primary beam is 87 / 95 arcsec 
at 33 / 30 GHz, corresponding to about 0.69 / 0.73 Mpc at 
the redshift of MRC 2104-242 / MRC 0943-242. 

High-2 radio galaxies are generally located in 
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Figure 3. a). 33 GHz radio continuum map of MRC 2104-242 (contour levels: 0.1, 0.3, 1.0, 1.7, 2.4 mjy bm"!). The cross indicates the 
location of the radio host galaxy, h). spectral line profile against the centre of the radio host galaxy. Shown is the full 1 MHz velocity 
resolution of CABB across the 2 GHz band. The x-axis shows the velocity in the rest- frame of the radio host galaxy (see Sect l2.lT3ll . 
The zoom-in shows a portion of the CABB data with the approximate bandwidth coverage of the old pre-CABB ATCA system (2 X 128 
MHz), c). rms noise per 1 MHz channel in the region of the radio source, d). Same as figure b, but data binned to 4 MHz channels 
(i.e. similar to the pre-CABB system; the zoom-in therefore gives a good representation of the data that could be obtained with the old 
2 X 128 MHz pre-CABB backend). e). rms noise per channel of 4 MHz in the region of the radio source. 

CABB (Manning smoothed to 2 MHz) 



-24°2B'30" 




Figure 4. a). 30 GHz radio continuum map of MRC 0943-242 (contour levels: 0.4, 1.0, 2.0, 3.0 mJy bm~^). The cross indicates the 
location of the radio host galaxy, h). Off-nuclear spectral line profile of redshifted CO. The spectrum is taken at the location marked by 
the arrow in figure aj and Hanning smoothed to a velocity resolution of 2 MHz across the 2 GHz CABB band. The x-axis shows the 
optical barycentric velocity in the rest-frame of the radio host galaxy (see Sect l2.1^ . c). rms noise per channel in the Hanning smoothed 
data of figure b) across the CABB band, derived across the central region. The arrow marks the rms noise level at the velocity of the 
tentative CO detection (right plot), d). Zoom- in of figure b, sho wing the tentative off-nuclear CO detection. The arrow indicates the 
redshift of the deep Lya absorption of H I gas jjarvis et al.|l2 003fl at which we centred our zero- velocity. The range of velocities of the 
emission-line gas in the giant Lya halo llVillar-Martfi7^tTul200a ) is also indicated in the plot. e). separate data-sets of the H75 and 
H168 array observations, both showing the tentative CO signal (for illustration purposes, the x-axis of the H75-array data is scaled-down 
by 3 mJy bm~^ in this plot). 



proto-cluster environ ments (e.g. IPentericci et al.l l2000al : 
IVenemans et all 120071 ). MRC 0943-242 is known to be lo- 
cated in a proto-cluster with many nearby compan ions de- 
tecte d in LyQ and with known redshifts ( Venemans et al.l 
l2007l ). There are 12 known LyQ companions within the 



primary beam and observing band of our obs ervations 
(Ivan Breukelen et al. I l2005l : IVenemans |2007| . see also 
Fig- El- None of these galaxies shows a clear CO detection 
above a 3a" limit, after correcting for primary beam attenu- 
ation. 
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Figure 6. Continuum image of MRC 0943-242. Siiown is an area 
the size of the FWHM of the ATCA's primary beam at 30 GHz. 
Flux densities are corrected for the primary beam response pat- 
tern, which causes the effective noise level to rise significantly 
towards the edge of the image. Contours: -0.4, -1.0, -2.0 (grey); 
0.4, 1.0 2.0, 3.0 (black/white) mjy bm~^. The red circles indi- 
cate the 12 Lya-bright galaxies in the field of MRC 0943-242 
(a nd with redshifts within th e CA BB observing band) de tected 
by Ivan Breukelen et al.l ||2005| ') and I Venemans et al] l l2007tl . 



Lya observations of the environment of MRC 2104-242 
are lacking and hence the cluster properties are unknown. 
No CO was detected within the primary beam above So" 
after correcting for primary beam attennuation. 



4 DISCUSSION 
4.1 H2 masses 

CO is an excellent tracer of molecular hydrogen, because the 
rotational transitions of CO are excited primarily by colli- 
sions with H2. A standard conversion factor Ox = Mh2/Lco 
[Mq (K km pc^)~^] is generally used to calculate the 
mass of the cold molecular gas (where Mh2 includes a 
fraction of the molecul ar gas that is in the form of he- 
lium - see for example ISolomon fc Vanden BoutI I2OO5I . for 
a review). For ultra -l umin ous infra-red galaxies (ULIRGs), 
iDownes fc SolomonI (|l998l ) derived a conversion factor of 
ax ~ 0.8 Mq (K km s~^ pc^)~^. This is in agreement 
with other observ ations of ULIRGs (|Solomon et al ] 1 19971 : 
lEvans et all I2OO2I ) as well as high-z sub-mm and star - 
forming galaxies (|Tacconi et al.l l2008l : IStark et all |2008| ). 
which imply that Ox ~ 0.8 — 1.6 M© (K km s~^ pc^)~^. 
We adopt a value of = 0.8 M© (K km s~^ pc^)~^ also for 
the two high-z radio galaxies that we study in this paper. 
We note, however, that there is a significant uncertainty 
in this conversion factor, since values as high as Qx ~ 5 
have been der i ved f o r molecular clouds in the M ilky Way 
i Scoville et"al] Il987l: IStrong et al.lll988l. see al s o iDickmanl 
1 19781) , iBloemen et all (|l986D . ISolomon et all (| 19871 )) as 



well as other nearby spi ral galaxies jPickman et al.l Il986l : 
ISolomon fc Barrettlll99ll ). 

Based on our 3cr upper limits on L'qq and assum- 
ing Qx = 0.8, we estimate that Mh2 < 2 x 10^" Mq for 
MRC 2104-242 and Mh2 < 6 x 10^° Mq for MRC 0943-242. 
The tentative off-nuclear CO detection in MRC 0943-242 
has an estimated molecular gas mass of Mh2 = 6x 10^° Mq. 

4.2 Molecular gas properties of high-z radio 
galaxies 

The upper H2 mass limits that we derive for MRC 2104- 
242 and MRC 0943-242 are comparable to H2 masses 
derived from CO de t ection s in high-g radio galax- 
ies (e.g. Scoville et al.l [199^: Papadopoulos et al] 200C ; 



De Breuck et al.[ 



■2003al lbl. 120051: iKlamer et al.1 120051 
Nesva dba et al. 20091. see also Isolomon fc Vanden Bout! 



2005. ); iMilev fc De Bre"uc5 (|2008l ) for reviews). However. 



as discussed in Sect. [T] most of these observations have 
targeted the higher rotational CO transitions, which could 
underestimate the total molecular gas content in these 
systems. CO (1-0) detections have been clai med for two 
high-2 radio galaxies, namely 4C 60.07 (2 = 3.8 Greve et akl 
I2OO4I ) and TNJ 0924-2201 (2 = 5.2 iKlamer et alj |2005| ). 
both with Mh2 = 1 x IO^'^Mq. Our derived upper limit 
on molecular gas mass in MRC 2104-242 (z = 2.491) is a 
factor 5 lower than this. 

Sub millimetre galaxies (SMGs) are likely merging sys- 
tems with a short-lived burst of extreme star formation and 
are believe d to be t he progenitors of local massive ellip- 
ticals (e.g. I Greve et~ al. 2005; Tac coni et al. 2008). In this 
sense, high-z radio galaxies and SMGs could be the same 
type of objects that differ on ly in their level of AGN activ- 
ity (e.g. iReuland et al.ll2007l ). although Unison et al.1 (|2008l ) 
argue that the violent AGN activity may occur predom- 
inant ly during th e early evolutionary stages of these sys- 
tems. Greve et al.l l|2005l ) derived a median cold gas mass of 
(M m) = 3.0x10^°M q among 12 SMGs detected in CO (see 
also iNeri et al.|[2003l ) . This is of the same order as the upper 
limits that we derive for the mass of cold gas in MRC 2104- 
242 and MRC 0943-242. 

Our derived upper limits on the molecular gas mass 
of MRC 2104-242 and MRC 0943-242 are lower than the 
H2 mass estimates for a non-negligible fraction of normal 
massive star forming gal axies at z ~ 1 — 2, d erived from 
CO (3-2) observations by iTacconi et all l|20ld . even when 
accounting for the much larger CO-to-H2 conversion fac- 
tor that they used). A similar result is seen by comparing 
the upper limits on CO in samples of h igh- 2 radio galaxies 
llEvans et al.lll996l; Ivan Oiik et aLlll997l ) with the results of 
iTacconi et all (|2010l ~ Confirmation by observations of larger 
samples in the same CO transitions might indicate impor- 
tant differences in molecular gas fraction, excitation proper- 
ties or chemical enrichment processes between high-2 radio 
galaxies and distant massive star forming galaxies. 

The H2 mass Umit of MRC 2104-242 is only a factor 
3 higher than the H2 content of the most CO-bright radio 
galaxies in the low redshift Universe, as studied from CO(l- 
0) observatons of a la rge sample of IR-bright radio galaxies 
bv lEvans et al. I (I2OO5I . corrected for the difference in the used 
Qx-value and cosmological parameters). The vast majority o f 
the low-2 radio galaxies in the sample of lEvans et aLl(|2005h . 
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however, contain signif icantly less molecular gas. T his was 
recently confirmed by lOcana Flaguer et al.l (2010) with a 
large sample of low-z radio galaxies not selected on IR- 
properties, for which they derive a median H2 mass of only 
Mh2 = 2.2 X 10** Mq. 

We note that many high-z CO detections to date are 
case-studies of galaxies that were pre-selected based on their 
properties at other wavelengths, such as a large sub- mm dust 
content or high infra-red (IR) luminosity. Both at low- and 
high-z there appears to be a relation between the far-IR 
(FIR) and CO luminosity in differe nt types of galaxies (see 
lEvans et al.l [2005l:[G reve et al]l2005l . and references therein). 
Such a relation would indicate that (radio) galaxies with a 
FIR luminosity in the range of ULIRGs (Lfir > 10^^ L©) 
contain a CO luminosity similar to the upper limit that we 
derive for MRC 2104-242 {L'co ~ few x lO" K km s"^ pc^). 
From Spitzcr observations of MRC 0943-242 at 24, 70 and 
160/im ([Seymour ct al. 2007), we estimate an upper limit on 
the total IR lumin osity of Ltir < 2x10 ^^Lb when using the 
approximation by iDale fc Heloul (120021). Follo wing the IR- 
CO re lation found bv lEvans et al.l l|2005l ') and iGreve et al.l 
(|2005h . this IR limit corresponds to an average CO lumi- 
nosity roughly a factor 2 lower than the L'qq upper limit 
that we derive for MRC 0943-242. The lack of detectable 
amounts of CO gas in MRC 0943-242 is therefore not un- 
usual based on its IR properties, but it shows that unbiased 
CO(l-O) observations of high-z radio galaxies are becoming 
feasible. 

Systematic searches for various CO transitions in un- 
biased samples of high-z (radio) galaxies are necessary to 
objectively investigate the overall content of cold molecular 
gas in the Early Universe. Our results show that system- 
atic and reliable searches for the ground-transition of CO in 
high-z (radio) galaxies are becoming feasible with existing 
broadband facilities that can target the 20-50 GHz regime, 
such as the ATCA and EVLA. 

4.2.1 CO in the vicinity of MRC 0943-242? 

In this Section we briefly discuss the possible nature of the 
tentative CO detection in the vicinity of MRC 0943-242, 
which needs to be confirmed before a more detailed analysis 
is deemed suitable. 

The tentative CO detection 60 kpc NE of the host 
galaxy) may be associated with a nearby companion galaxy, 
althou gh no companion has been detected in Lya at tha t lo- 
cation ( van Breukelen et al.ll200"5l : IVenemans et al]|2007l ). so 



any such galaxy would have to be Lya-faint. Alternatively, 
the tentative CO detection may repres ent cold gas in the 
outer part of the qu i escen t Lya halo (|Villar-Martm et al.l 
I2OO3I '). iBinette et al.l (|2000l ') show that C IV absorption is 
associated with the deep Lya absorption in MRC 0943-242 
and derive that this reservoir of absorbing gas is also lo- 
cated in the outer halo (i.e. ouside the radio cocoon). If 
confirmed, the cold gas properties of MRC 0943-242 resem- 
ble thos e found in the high-z radio galaxies TXS 0828+193 
(z = 2.6lNesvadba et al.ll2009l ) and B3 J2330+3927 (z = 3.1 
iDe Breuck et al.ll2003al ). 

The only two known high-z radio galaxies in which 
CO(l-O) has been detected (4C 60.07 and TNJ 0924-2201; 
see Sect. 14. 2|) also show indications that the CO gas 
may not be aligned with the central location of the host 



galaxy (|Klamer et al.ll2004l : llvison et al.ll2008l '). In particu- 
lar 4C 60.07 shows an apparent deficit of molecular gas in 
the radio host galaxy, while CO appears to be present in a 
merging companion and associated tidal debris (jlvison et al.l 
,2008 ). If confirmed, a more detailed comparison between the 
CO(l-O) properties of these systems deserves further atten- 
tion. 

The position angle of the radio source in MRC 0943-242 
(which has a total linear size of about 4" an d is therefore 
unres olved in our observations) is PA = -74° (jCarilli et al.l 
[l993). This is roughly within 45° of the location of the ten- 
tative CO detection from the central region of the radio 
host galaxy. If confirmed, this may resemble alignments that 
[ki amer et all (I200J) argue exist among other high-z radio 
galaxies. 

4.3 Radio continuum 

Both radio sources are clearly detected in our sensitive 
(a ~ 30 /iJy) 7mm continuum observations. Their spectral 
indices are relatively steep from ~30 GHz down to 1.4 GHz, 
with no evidence for spectral curvature within this large 
range of frequencies. This indicates that there is no turn-over 
due to synchrotron losses or inverse Compton cooling up to 
~115 GHz in the restframe of these radio sources. This is 
consistent with continuum observations of a l arge sample of 
high-z ultra-steep spectrum radio galaxies bv lKlamer et al.l 
(|2006h . who also find relatively steep power law spectral 
energy distributions (SEDs) with no evidence for spectral 
steepening up to several tens of GHz in the rest frame. A 
detailed analysis of this phenomenon is crucial for under- 
standing the electron acceleration mechanism or environ- 
mental properties of high-z radio sources, but is beyond the 
scope of this paper. 



5 CONCLUSIONS 

We presented the first 7mm observations of two high-z ra- 
dio galaxies (MRC 2104-242 and MRC 0943-242) with the 
2x2 GHz Compact Array Broadband Backend. Our re- 
sults demonstrate the feasibility of using ATCA/CABB for 
spectral- line work at high redshift. We also presented 7mm 
continuum images of the two high-z radio galaxies, with a 
typical rms noise level of ~30 /iJy beam~^. The enhanced 
spectral-line and continuum capabilities of ATCA/CABB in 
the millimetre regime complement those of other large exist- 
ing and upcoming observatories, such as Pdbl, EVLA and 
ALMA. 

From our CO (1-0) data we derive upper limits on the 
H2 mass of Mh2 < 2 x 1O^°M0 for MRC 2104-242 and 
Mh2 < 6 x 10^°Mq for MRC 0943-242 (ax = 0.8). These 
upper limits are of the same order as H2 mass estimates de- 
rived from CO detections of other high-z radio galaxies and 
SMGs, but lower than the mass of molecular gas detected 
in a non- negligible fraction of normal star forming galaxies 
at z ~ 1 - 2. For MRC 0943-242 we also find a tentative 
CO (1-0) detection at about 60 kpc distance from the cen- 
tral region of the host galaxy, but this needs to be confirmed 
with additional observations. 

The spectral index of both MRC 2104-242 and 
MRC 0943-242 is relatively steep with a ^ -1.5 between 
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1.4 and 30 GHz. There is no evidence for spectral curvature 
up to ~ 115 GHz in the rest frame of these radio sources. 
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